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ABSTRACT 

Aims. The detection of small-amplitude nonradial modes in high-amplitude 6 Set (HADS) variables has been very elusive until at 
least five of them were detected in the light curve of V974 Oph obtained from ground-based observations. The combination of radial 
and nonradial modes has a high asteroseismic potential, thanks to the strong constraints we can put in the modelling. The continuous 
monitoring of ASAS 192647-0030.0eeCoRoT 101155310 (P=0.1258 d, 7=13.4) ensured from space by the CoRoT (Convection, 
Rotation and planetary Transits) mission constitutes a unique opportunity to exploit such potential. 

Methods. The 22270 CoRoT measurements were performed in the chromatic mode. They span 152 d and cover 1208 consecutive 
cycles. After the correction for one jump and the long-term drift, the level of the noise turned out to be 29 /jmag. The phase shifts 
and amplitude ratios of the coloured CoRoT data, the HARPS spectra, and the period-luminosity relation were used to determine a 
self-consistent physical model. In turn, it allowed us to model the oscillation spectrum, also giving feedback on the internal structure 
of the star. 

Results. In addition to the fundamental radial mode f\ =7.949 d~' with harmonics up to 10/i, we detected 12 independent terms. 
Linear combinations were also found and the light curve was solved by means of 61 frequencies (smallest amplitude 0.10 mmag). 
The newest result is the detection of a periodic modulation of the f\ mode (triplets at ±0.193 d -1 centred on f\ and 2/i), discussed 
as a rotational effect or as an extension of the Blazhko effect to HADS stars. The physical model suggests that CoRoT 101 1553 10 is 
an evolved star, with a slight subsolar metallic abundance, close to the terminal age main sequence. All the 12 additional terms are 
identified with mixed modes in the predicted overstable region. 

Key words. Stars: variables: 6 Set - stars: oscillations - stars: interiors - stars: individual: CoRoT 101 155310 

1. Introduction logical separation is based on an amplitude limit of 0.20-0.30 

mag, but actually it has never been considered as physically sig- 
High-amplitude 5 Set (HADS) stars occupy the central part nmcant . Moreover, HADS stars could not be further considered 
of the lower part of the instability strip between the horizontal mono _ or do ub le-m ode radial pulsators after detecting addi- 
branch and the main sequence (|McNamara||2000|). In the same tional modes Jp O rettill2003L and references therein). Their iden- 
region of the instability strip, low-amplitude 6 Set stars span a tification as nonrad i a l modes opened the possibility of consider- 
wider range of temperatures and luminosities. The phenomeno- ing HADS stars as new tools fo r asteroseismic studies. Ind eed, 
recent theoretical investigations ( Suarez et al.ll2006al l2007l) in- 

* The CoRoT space mission was developed and is operated by the die ate that rotation effects may provide useful additional infor- 

French space agency CNES, with the participation of ESAs RSSD and mation on the mode identification, in particular for the nonra- 

Science Programmes, Austria, Belgium, Brazil, Germany, and Spain. dia i modes ro tationally coupled with the fundamental and first- 

This work uses ground-based spectroscopic observations made with overtone radial modes That large amplitude modes could be 

the HARPS instrument at the 3.6m-ESO telescope (La Silla, Chile) i deri tifi ed radial allo wed us to use HA DS stars as distance indi- 

under the ESO Large Programme LP182.D-0356 and complementary ... i l hnno LFi > i I hnm k tu „■ 

. . t . " ~ , . „. . .~ ;. Cf t . cators Poretti et al. 2008; Gars et al. 2010 . The most luminous 

photometric measurements made at the Piszkesteto Mountain Station , TAI ^„ : — , — ', — , — , — : — — TT 

of Konkoly Observatory (Hungary) HADS stars are just below the horizontal branch, and the debate 

** Tables 1 and 2 are available in electronic form at the CDS whether the light curves of HADS stars could show similarities 



anonymous 



ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via with those of RR Lyr stars is still open (|Breger||2010|). 



http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/vol/page 
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A great improvement in the observational techniques and 
data quality was needed to investigate both the pulsational con- 
tent and the behaviour of the predominant radial modes. This 
opportunity could only be offered by a space mission. 

2. The CoRoT data 

The variability of the 13 th -mag star located at a — 
19 h 26 m 47.25, 6 = -00° 30' 01744 (J2000) was disc overed dur- 
ing th e ASAS survey (All Sky Automated Survey; iPoimanskH 
l2Q00h . The star is named 192647-0030.0 and classified as 
a 6 Set variable with P=0.1258 d in the ASAS catalogue 
(http://www.astrouw.edu.pl/asas/). This part of sky, 
i.e., the equatorial plane not far from the galactic centre, 
is well within one of the "eyes" of the CoRoT space mis- 
sion (COnvection, Rotation and planetary Transits; iBaglin et"ai1 
120061: lAuvergne et aT1l2009l) . The inclusion of ASAS 192647- 
0030.0 in the exoplanetary field of the first long run in 
the direction of the galactic centre (LRcOl) was done afte r 
searching the catalogue compiled by iBenko & Csubrvl (120071) . 
Therefore, we specifically requested the CoRoT monitoring of 
ASAS 19264 7-0030.0 as part of the Additional Programmes 
(IWeissI 120061) . In the CoRoT database the star was named 
CoRoT 101 1553 10sUSN OA2 0825-15821309 (V=13.43, B - 
V=+0.60, V - i?=0.21; iDeleuil et al.1 |2009h . There is no 
relevant contamination from nearby stars, since the bright- 
est star included in the CoRoT mask is 2.9 mag f ainter 
than CoRoT 101155310 in V light dDeleuil et al.l 12009b . The 
"CoRoT Variability Cla ssifier" automated supervised method 
dDebosscher et al.l l2009t) classified it as an RR Lyr star only 
on the basis of the Fourier parameters. Since some ov erlap is 
prese nt with the HADS stars in this parameter space dPorettil 
1200 lb . the classification system has been consequently revised 
to better separate RR Lyr from HADS stars. 

The LRcOl started on May 15, 2007 and finished 152 d 
later, on October 14, 2007. This means that CoRoT 101155310 
was monitored for 1208 consecutive cycles. The exposure time 
in the CoRoT exoplanetary channel was 512 s. This value re- 
mained constant over the long run, since no regular feature sug- 
gesting planetary transit was noticed. The Heliocentric Julian 
Dates (HJD) were reported to the times at half exposure. The 
very high duty cycle allowed us to collect 25597 datapoints. 
CoRoT 101155310 was observed in the chromatic mode, by 
means of a bi-prism located in the focal block of the instru- 
ment, b etween the dioptr ic objective and the plane of the de- 
tectors dBarge et al.ll2006l) . The image is a small, low-resolution 
spectrum whose area is about 300 arcsec 2 on the sky. The pho- 
tometry is performed in three different colours, roughly red (rc), 
green (gc), and blue {be), to supply a tool able to distinguish 
between planetary transits (colour independent) and stellar ac- 
tivity (changing with the wavelength). The white light, i.e., the 
sum of the three contributions, is more suitable for signal detec- 
tion thanks to its larger number of counts. CoRoT 101155310 
was also observed with the 1.0 m telescope of the Konkoly 
Observatory at the the Piszkesteto Mountain Station (Hungary) 
on two nights in June 2009. CCD photometry with standard BY I 
filters was planned to compare the amplitudes in the CoRoT 
passbands with those in standard ones. The 88 measurements 
in each filter cover one complete pulsational cycle, and we could 
also determine one time of maximum brightness. 

The CoRoT measurements are a sort of absolute photometry, 
since the measured flux is not corrected with respect to any ref- 
erence source or star. Therefore, instrumental changes and drifts 
have a direct influence on it, on both short- and long-term scales. 



Different detrending algorithms can be used, based on moving 
means or polynomial fits. After several trials, we were more sat- 
isfied by an approach based on a preliminary frequency analy- 
sis. Indeed, the main frequency f\ =7.949 d _1 and other high- 
amplitude terms were immediately detected in the raw data. We 
calculated the least-squares fit of the raw data by using a prelimi- 
nary solution based on the terms with the largest amplitudes. The 
residuals in white light (Fig. Q] upper panel) show a long-term 
drift and more rapid and apparently irregular changes around 
JD 2454253 and 2454267. Moreover, in the second half of the 
observations the distribution of the outliers strongly changed 
and many points appeared well below the others. This so-called 
"rain" of points from the middle strip has an instrumental ori- 
gin because it is observed in other stars belonging to the same 
pointing. To correct the irregular behaviour of the mean mag- 
nitude we calculated the least-squares fits on four consecutive 
cycles. These fits returned the 302 values of the mean magni- 
tudes (Fig.Q] upper curve in the middle panel), which mimics the 
trend of the residuals. The standard deviations of the 302 solu- 
tions range from 0.004 to 0.016 mag. To remove the instrumental 
variations, we interpolated the values of the mean magnitudes at 
the time of each individual measurement of CoRoT 101155310 
and then subtracted it. 

After that, we recalculated the least-squares fits on the 302 
bins and obtained very homogeneous values of the mean magni- 
tudes (Fig.Q] lower curve in the middle panel). By applying this 
procedure we also removed any oscillation around 0.5 d, but as 
we see later (Sect. 13.4b this did not have a strong impact on our 
solution. At this stage we also removed the outliers. A post-facto 
analysis of the raw and the cleaned datasets showed that the fre- 
quency content of the two timeseries is the same. The bottom 
panel of Fig.Q]shows an example of the excellent quality of the 
final timeseries and also allows us to appeciate the continuous 
CoRoT monitoring. 

The same technique was applied to the coloured data. The 
trend in rc light is roughly the same as in white light (top panel 
of Fig.Q}, since the total flux in this light (about 182000 adu) is 
much greater than of gc (about 22600 adu) and be (about 32300 
adu). The main difference is in the decline after JD 2454310, 
which is sharper in rc light than in white light. Actually, the 
fluxes observed in the gc and be light increased after this date, 
and what we measured in white light is the total flux variation. 
This changeover concerns those days when the "rain" effect ap- 
peared, and it is very probably a jitter effect through the photo- 
metric mask of CoRoT 101155310. 

3. The frequency content 

The frequency analyses were performed on different versions of 
the CoRoT timeseries (raw data, detrended data, coloured data, 
with and without outliers, etc.) and they supplied very simi- 
lar results. The main change is in the level of the noise and 
consequently in the threshold for the significance of the peaks. 
Therefore, for the sake of clarity, we discuss here the detection 
performed in white light. The frequency values were firstly ob- 
tained b y means of the iterative sine-wave, least-squares fitting 
method (Vani £eklll971h an d then refined by the MTRAP algo- 
rithm (Carpin g et al.lll98~7l) . both at intermediate stages and at 
the end of the whole analysis. 

The level of the noise in the final timeseries (22270 data- 
points in white light) is 29/imag in the 0-100 d -1 region. The 
long-term drift left a slightly higher noise in the 0-5 d -1 region, 
namely 41 //mag. We accepted 61 peaks as having a stellar 
origin, i.e., not counting the orbital frequency of the satellite 
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Fig. 1. Top panel: curve of the residuals (a preliminary solution was 
subtracted from the original data) showing the long- and short-term 
instrumental effects. Middle panel: behaviour of the mean magnitudes 
before (upper curve) and after (lower curve) the removal of the instru- 
mental effects shown in the upper panel. Each bin spans four pulsational 
cycles. Bottom panel: the example extracted from the final timeseries 
of CoRoT 101155310 shows the continuous monitoring obtained from 
space by CoRoT. 



/or*=13.97 cT , the sidereal frequency /, ( /=1.0027 cT , and their 
harmonics and combinations. To do that, we firstly considered 
the peaks with an amplitude greater than 0.12 mmag, i.e., with 
a signal-to-noise ratio (SNR) above 4.0. Then, we noted that 
several combination terms (harmonics and coupling terms) were 
also detected below this limit. Therefore, we relaxed the min- 
imum amplitude to 0.10 mmag, i.e., 3.5 times the noise level. 
Some combination terms (e.g., /1+/2+/3, 4/1-/4, 3/+/) are 
just below the acceptance threshold. 

Figure |2] shows the steps in the solution of the CoRoT light 
curve of CoRoT 101155310, which led to the detection of four 
categories of frequencies: 

1. the predominant term / and its harmonics; 

2. the independent terms, from / to /13; 

3. the modulation term f,„ of the term/; 

4. the combination terms. 

The parameters of the least-squares solution (cosine series) 
based on the 61 frequencies thus detected are listed in Table Q] 
The residual rms is 0.0027 mag. The values of the har- 
monics and combination terms were obtained from the val- 
ues of the independent frequencies and from the modula- 
tion frequency /m=0.1933 d' 1 , fo r which formal errorbars 
dMontgomerv & 0'Donoghiuelll999l) are given between brack- 
ets. This "locked" solution was calcu lated by means of the 
MTRAP algorithm dCarpino et al.lll987l) . The residual spectrum 
(Fig. |2] panel e) shows some peaks at the values of the linear 
combinations of f or b— 13.97 d _1 and /,</= 1.0027 d~'. These spu- 
rious peaks are due to the passage of the satellite over the South 
Atlantic Anomaly, which occurred twice daily. The environmen- 
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Fig. 2. Subsequent steps in the detection of frequencies in the ampli- 
tude spectra of CoRoT 101155310. 



tal conditions (e.g., the effects of the eclipses on the electronics 
units, the eclipse duratio ns, the difference in th e Earth's albedo 
of the overflown regions: IXuvergne et al. 2009) are affecting the 
CoRoT photometry in s uch a complicated w ay that, as in the case 
of CoRoT 101128793 (iPoretti et al.ll2010t) . the usual technique 
of prewhitening was not able to clean the data in a satisfactory 
way. 

3. 1 . The / term and its harmonics 

The peaks at / =7.949 d -1 , 2/, and 3/ clearly stand out in the 
amplitude spectrum of the original data (Fig. |2 panel a). The 
spectral window is free of any relevant alias structure, the high- 
est peak being at the power level of 0.005 (see the inserted box). 
We detected harmonics up to 10/ in the CoRoT light curve, i.e., 
about twice those detected in ground-based timeseries. Since the 
amplitude of harmonics is monotonically decreasing, the light 
curve is smoothed (Fig. [3] upper panel), with out the bump s ob- 
served in RR Lyr stars, both monoperiodic (lPorettill200ll) and 
multiperiodic (IChadid et aIll2010tForetti et al.boiOl) . 

The Fourier parameters (021 = 4.19 rad, ifoi = 2.04 
rad, 04i = -0.27 rad, and ^21 = 0.38) are in excellent 
agreement with those calculat ed from the OGLE light curves 
of HADS stars dPorettil l200lh . These parameters suggest that 
CoRoT 101 155310 is pulsating in the fundamental radial mode. 

We calculated the times of the maximum brightness by 
fitting the measurements of the well-covered cycles with 
a cubic spline. In such a way we could determine 1090 
times of maxima (Table 2) and in turn the ephemeris 
Max = HJD 2454236.7014 +0.12579931 • E 

+0.0001 +0.00000011 
by means of a least-squares line. The O-C values (differences 
between the observed and calculated times of maxima, also 
listed in Table 2) were determined by using this ephemeris. 
The CCD BVI measurements performed at the Piszkesteto 
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Table 1. Identification, frequencies, Fourier amplitudes, and phases (To=HJD 2454236.7502) of the terms detected in the data of 
CoRoT 101155310 (22270 datapoints). 



ID 


Frequency 


Amplitude 


Phase 


ID 


Frequency 


Amplitude 


Phase 


m 




[d- 1 ] 


[mag] 


[0,2*] 




[d- 1 ] 


[mag] 


[0,2*] 







Main frequency and harmonics 






Additional terms identified as independent modes 






7.949172(1) 


0.137632(26) 


5.3775(2) 


A 


10.847014(015) 


0.006224(26) 


6.0785(0041) 


0.733 


v\ 


15.898343 


0.041859 


2.3933 


A 


10.479540(016) 


0.005689(26) 


5.8511(0045) 


0.759 


3/. 


23.847515 


0.014404 


5.6197 


A 


10.690302(030) 


0.003076(26) 


5.8388(0083) 


0.744 


4/, 


31.796687 


0.004433 


2.4116 


A 


9.532648(075) 


0.001230(26) 


3.1547(0208) 


0.834 


5/, 


39.745859 


0.002146 


5.4146 


A 


13.175878(147) 


0.000629(26) 


1.1226(0407) 


0.603 


6/. 


47.695030 


0.001150 


2.2493 


A 


9.796144(340) 


0.000273(26) 


5.7221(0937) 


0.811 


7./, 


55.644203 


0.000552 


5.3327 


A 


10.203291(385) 


0.000241(26) 


4.1313(1062) 


0.779 


8,/, 


63.593373 


0.000317 


2.3551 


A 


8.232572(450) 


0.000206(26) 


0.2385(1242) 


0.966 


9/. 


71.542544 


0.000130 


5.3127 


/io 


12.635143(455) 


0.000204(26) 


0.0702(1254) 


0.629 


10/, 


79.491717 


0.000100 


2.6673 


/i, 


13.224629(504) 


0.000184(26) 


3.3973(1391) 


0.601 










Ai 


9.325259(549) 


0.000169(26) 


5.2843(1514) 


0.852 










fa 


8.354265(672) 


0.000138(26) 


1.0701(1854) 


0.952 



Modulation terms with f„,=0. 1933 d 1 Linear combinations between / and/2 



f-fm 


7.755895 


0.000194 


0.1805 


A +A 


18.796186 


0.002261 


2.2225 


f\ + fm 


8.142448 


0.000145 


1.9766 


A-A 


2.897843 


0.001524 


3.6718 


2/l ~fm 


15.705067 


0.000131 


2.5886 


2/1 +A 


26.745358 


0.001323 


5.5501 


2/l + f m 


16.091183 


0.000106 


4.5329 


2/1-/2 


5.051329 


0.000814 


5.5587 










3/i +A 


34.694529 


0.000687 


2.4496 


Linear combinations between / and A 




3/,-/2 


13.000501 


0.000219 


2.8031 










4/i +A 


42.643701 


0.000344 


5.5241 


/ + A 


18.428712 


0.002247 


2.1349 


5/ +A 


50.592873 


0.000213 


2.6688 


A- A 


2.530368 


0.001398 


3.4160 


5 A -A 


28.898844 


0.000131 


4.2861 


2/1 +A 


26.377884 


0.001 104 


5.3307 


6/1 +A 


58.542046 


0.000144 


5.5907 


2/1-/3 


5.418803 


0.000781 


5.8036 










3/1 +A 


34.327056 


0.000615 


2.2546 




Linear combinations between / and / 5 




3/1-/3 


13.367975 


0.000169 


2.7857 










4/1 +A 


42.276226 


0.000303 


5.3800 


A +A 


17.481819 


0.000322 


5.6402 


5/1 +A 


50.225399 


0.000199 


2.2301 


A-A 


1.583476 


0.000145 


0.5953 










2/. +A 


25.430991 


0.000139 


2.4076 


Linear combinations between / and / 4 






Linear combinations between / and f 6 




2/4 


21.380604 


0.000164 


2.8776 


A +A 


21.125050 


0.000217 


3.6347 


/1+/4 


18.639474 


0.001137 


2.0182 


A-A 


5.226706 


0.000141 


4.2040 


A-A 


2.741130 


0.000813 


3.4594 


2/1 +A 


29.074221 


0.000153 


0.1586 


2/ +/4 


26.588646 


0.000620 


5.3297 


2/1-/6 


2.722466 


0.000135 


5.2315 


2/1-/4 


5.208042 


0.000479 


5.6289 










3/i +A 


34.537817 


0.000309 


2.1560 




Linear combinations between ft and / 




3/1-/4 


13.157213 


0.000127 


2.6059 










4/ +/4 


42.486989 


0.000121 


5.4129 


A+A 


21.326554 


0.000194 


3.3195 


5/i +A 


50.436162 


0.000105 


2.3231 











Observatory in 2009 provide an excellent validation. The max- 
imum in / light was observed at HJD 2455011.4997 when 
the above ephemeris predicts HJD 2455011.4994+0.0007 for 
E=6159. The light curve of this maximum is flatter than those 
of the B and V curves, whose maxima were observed at 
HJD 2455011.4933 and HJD 245501 1.4910, respectively. 

3.2. The independent terms 

Twelve independent frequencies were found in the interval 8.23- 
13.22 d _1 , i.e., at periods shorter than that of the fundamental 
radial mode f\ (Table [TJ. These additional modes could not be 
detected in the ASAS timeseries since the level of the noise is 
much higher (0.020 and 0.025 mag in V and /, respectively) than 
their amplitudes (0.006 mag or less). 



A close, but non-equidistant triplet of frequencies 
(/2=10.48, / 3 =10.69, and / 4 =10.85 d" 1 , Fig. |2] panel b) 
appears at an amplitude level similar to that of Af\, i.e., in the 
3.0-6.2 mmag range. Is CoRoT 101155310 a double-mode 
HADS star? The frequency ratio between the fundamental 
and first overtone modes in HADS stars is well constrained 
in the 0.771-0.7 74 interval (see Table 2 and Fig. 4 in 
iPoretti et al.1120051) . In turn, this implies that the first overtone of 
CoRoT 101155310 is expected in the interval 10.27-10.31 d -1 . 
None of the terms of the triplet matches this value and so the 
triplet is composed of nonradial modes. 

The amplitude of f\ is half that of fa, and this ratio is also 
maintained for fs, f(, (Fig. [2] panel b), and fi. After fi (am- 
plitude 0.27 mmag), the amplitude of the subsequent terms is 
slightly decreasing and we detected four other modes (fi,fg, fio, 
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Fig. 3. Top panel: the CoRoT data folded with the main pulsational 
frequency f\ (upper curve) and the residuals after prewhitening of the 
61 frequencies (lower curve). Bottom panel: amplitudes of the f\ term 
calculated from measurements spanning four consecutive pulsational 
cycles. 



and f\\) in less than 0.10 mmag (Fig. [2 panel c). Still, no reg- 
ular separation is observed among these peaks. More generally, 
we did not find the same separation between two terms from f\ 
to /13. Regarding the identification of the first radial overtone, 
the fq term provides the ratio closest to the expected one, but it 
appears too discrepant (see also Sect. [6). 



3.3. The modulation of the fi component 

The side peaks of /1 (7.7559 and 8.1424 d" 1 ; Fig.El panel d) and 
2/i (15.7051 and 16.0912 d~') are the most surprising features 
in the amplitude spectrum of CoRoT 101 155310. They form an 
equidistant triplet with a frequency separation /,„=0.193 d _1 . In 
practice, this triplet is the signature of a modulation of the main 
component f\ . The amplitudes of the modulation terms are much 
smaller than those of the independent terms discussed above. 
This hampers any direct detection of the periodicity in the O-C 
or the magnitude-at-maximum diagrams, since the effects of the 
independent modes mask those of the modulation, also taking 
the long (512 s) CoRoT exposures into account. The residual 
signal left around f\ (noted as "Long term" in Fig. [2] panel d) 
will be discussed in Sect. 17.21 

The two frequencies /i2=9.325 and /o=8.354 d are 
close, but not coincident, to the combination terms f$- 
/ ffl =9.340 d _1 and /i+2/„=8.335 d _1 . The hypothesis that they 
are actually combination terms implies an f,„ value significantly 
different from 0.193 d" 1 (0.208 and 0.202 d -1 , respectively). 
Therefore, we should admit that the modulation is not the same 
or that these two frequencies are independent modes. Since the 
occurrence of a combination term such as fs-f m is quite unusual 
we prefer the latter explanation. 



-30.0 -20.0 -10.0 0.0 10.0 20.0 

b _ r -$ r fdegl 

Fig. 4. Phase shifts and amplitude ratios of the frequencies from fi to 
as obtained from CoRoT coloured timeseries. 



3.4. The combination terms 

The amplitude spectrum of CoRoT 101155310 is also rich in 
combination terms (Fig. [2] panel b), as expected from the am- 
plitude of the independent terms from f\ to fs, greater than 
1 mmag. The terms involving differences were also detected, and 
the lowest frequency is /5-/i=1.58 d _1 . We can argue that the 
bin on four pulsational cycles (i.e., 0.5 d) did not hamper detec- 
tion of peaks below 2 d -1 , if any. Patterns of combination terms 
with the harmonics up to 5/i are detected between / and f2, f\ 
and/,/ and /4, up to 2/i with/5 and/,. Only four frequencies 
(3/i-/2=13.00 d~\ 3/i-/ 3 =13.37 d" 1 , 3/-/t=13.16 d- 1 , and 
/i+/ 5 =17.48 d~') are within the interval 7-18 d _1 . We consid- 
ered them as combination terms rather than independent modes, 
since they follow the patterns described above closely (see also 
Table[T]i. We also detected the term fi+fi, i.e., the one involving 
the additional terms with the largest amplitudes. The modulation 
does not affect the combination terms; i.e., we did not detect any 
term having the form / +/„,+/ . 



4. The fit of the light curves in the different CoRoT 
passbands 

Calculation of the phase differences and amplitude rati os in dif- 
feren t passbands could allow the mode identification (iGarridol 
120001) . Unfortunately, the CoRoT photometric system is not a 
well-defined one , depending on the star itself (see Sect. 7 in 
iBarge et aT]|2006l) . To perform a comparative procedure, we cal- 
culated the least-squares fits of the coloured CoRoT and BVI 
Konkoly timeseries. The full amplitudes of the light variations of 
CoRoT 101155310 are 0.432, 0.385, and 0.273 mag in the b c , 
gc, and rc CoRoT passbands and 0.459, 0.367, and 0.215 mag in 
the BVI filters, respectively. Therefore, the be and gc passbands 
are not separated well since both are between the B and V ones, 
while the be and rc passbands are both a little redshifted with 
respect to the B and V ones. We plotted the phase shift <pi,- r - 4> r 
vs. the amplitude ratio Ab- r /A r . The position of the mode / is 
very far from those of the modes fj—fe, all clustering in a narrow 
region (Fig. |4j. Positive values of the phase differences suggest 
( — (i.e., radial) modes, negativ e values are likely I > 1 (i.e., 
nonradial) modes dWatsonlll988l) . Moreover, we can see very 
good agreement between the observed positions of the modes 
h ~f(> and the th eoretical locus of the 1-2 modes (see Fig. 18 
in IWatsonl[l988l) . The identification of these frequencies with 
the 1-2 modes must be considered as only preliminary, due to 
the uncertainties in the characteristics of the CoRoT photomet- 
ric system. 
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The phase errorbars associated to are already too large 
to clearly point out the mode identification, and this becomes 
impossible for the frequencies from fi to /13. In particular, we 
cannot obtain clues on the identification of fg as the first radial 
overtone. 



5. The stellar physical parameters 

The spectrum of CoRoT 101155310 was obtained with the 
HARPS instrument mounted at the 3.6-m ESO telescope on the 
nights of 29-30 June 2009. Due to the star's faintness and the 
short periods, the resulting SNR value was around 30 (/?=80000, 
exposure time 2300 s), and we could use the spectrum to esti- 
mate the physical parameters only. 

The mean line p rofile was computed with the LSD software 
(iDonati et al] 1 19971) and the Fourier transform of such profile 
supplied y e(? sin /=26+ 1 kms . However, this determination is 
affected by the broadening owing to the radial pulsation of the 
fx mode. Since we do not know the amplitude of the radial ve- 
locity curve of CoRoT 101 155310, we estimated it by using the 
2KI k. V value. For radial pulsation it is ar ound 100 kms 1 mag 1 



Table 3. Observed and physical parameters of 
CoRoT 101 155310 determined in this paper. 



(ISmithlll982l: iPoretti & Antonelloll 19881). the highest valu e was 
probably measured for V356 Aur ( Poretti et al.l 11987k i.e., 



of 



130 kms mag . Since AV = 0.37 mag in the case 
CoRoT 101155310, we got 2^=37 kms 1 and 48 kms 1 , re- 
spectively. We also estimated the 2K value in the specific case 
of CoRoT 101 155310, obtai ning 2K=12 km s' 1 using the fheo - 
retical approach described bv lMova & Rodrfguez-L6pezl(l2010l) . 
The ephemeris determined in Sect. 13.11 allowed us to establish 
that the HARPS spectrum started after the maximum light and 
ended around the halfway point of the descending branch. In 
this interval the radial velocity due to pulsation is getting close 
to zero after its maximum value. We evaluated the effects of 
both the radial pulsation and the exposure time on the mean lin e 
profile by using the software package FAMIAS (IZimall2006l) . 
Our simulations indicated that the 2K values determined above 
could produce a maximum broadening of 5 kms -1 in the de- 
termination of the V eq sin i value. By combining all these re- 
sults, we accepted 23+4 kms -1 as the value of the V^sim 
ofCoRoT 101155310. 

We computed a grid of synthetic spectra with the /?=80000 
resolution and the observed line broadening by means of 
the A TLAS9 and SYNTHE codes dKuruczlll993"l: iHeiter et alj 
120021) . The grid covers the ranges 7000< T eff <8000 K, 3.2< 
logg <4.8 , and -0.5<[Fe/H |<1.0. A microturbulent velocity of 
2 kms 1 dGrav et ail 1200 lb was considered. We compared the 
observed spectra with the synthetic grid in chosen wavelength 
regions in order to estimate the physical parameters of the star: 
the Ha and HB wings for r eff and the Mg triplet (5167, 5173, 
and 5184 A) for logg. We obtained finally T eS = 7300;™ K 
and logg = 3.75 ± 0.20 dex. These values were t hen used for 
a com plete abundance analysis (see Sect. 4.2 in IPoretti et alJ 
120091) . which revealed a slight metallic solar un derabundance, 
i.e., |E e/H1=-0.15±0.20 dex ([Fe/H] Q =7.54 dex; lAsplund"etai1 
120091) . 



Identifying the fundamental radial mode allows us to cal- 
culate the absolute magnitude. We get M y-1.46+0.10 by us - 
ing the period-luminosity relation given bv iPoretti et ai1(l2008l) . 
Other relations give similar results for log P = -0.900 (see 
Fig. 8, ibidem). The values My-=1.46 and T eff = 7300K puts 
CoRoT 1011 55310 in the mid dle of the HADS instability strip 
(see Fig. 2 in lMcNa mara 200(|. 
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13.94 


logg 
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+0.20 


[Fe/H] c 


-0.15+0.20 dex 


K* 


12.11 


V eq sin; 


23+4 kms" 1 


(V - K) 


+0.70 


L 


19.8+1.8 L e 


-Pfimd 


0.1258 d 


R 


2.8+0.4 R Q 


My 


1.46±0.10 


M d 


1.71+0.02 Mq 


M b0] 


1.49+0.10 
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1485+10 Myr 



a From lSkrutskie et al.l J20061) 

b Pulsational models suggest r eff =7000±50 K 

c Pulsational models suggest [Fe/H] =-0.30+0. 20 dex 

d Values from pulsational models 



The approach described by iBreger & Bregmanl (Il975l) sup- 
plied the relation 



log Q = -6.456 + logP + 0.5 logg + 0.1M bol + log r eff . 



(1) 



We know that P=0.126 d from the CoRoT observations and that 
2=0.033 d from the identification of P as the fundamental ra- 
dial mode. Sinc e the bolome tric correction is +0.03 mag for 
T eS = 7300K ( lTorresll2010l) . we obtain M bo i=1.49+0.10 and 
L = 19.8 ± 1.8Lq from the M V -IA6 value calculated above. 
The unknowns g and r e ff must therefore satisfy the condition 



logg + 2.01ogr eff = 11.45(+0.01). 



(2) 



This condition holds for the central values of log g and T e g ob- 
tained from HARPS spectra and proves the self-consistency 
of the physical model. We can then calculate the radius R = 
2.8 + 0.4 R Q and the mass M = 1.6 + 0.7 Mq from the L, T eff , 
and \ogg values. 

We further checked the consistency of the physical model 
by combining our BV photometry with the 2MASS photom- 
etry. The mean BV magnitudes are < V >=13.37 and < 
B >=13.94. These values are in good agreement with those 
provided by the EX ODAT catalogue (V=13.43, B - V=+0.60; 
iDeleuil et al1l2009l) . Indeed, the epoch of the EXODAT obser- 
vations corresponds to the phase 0.23 p, i.e., roughly the mid- 
dle of the descending branch (Fig. |3). We calculated the equiv- 
alent widths (EWs) of the Nal lines at 5890.0 and 5895.9 A 
from the HARPS spectrum. We could estimate the interstellar 
redd ening Eb-v=+Q-2Q by m eans of the Eb-v~EW relation- 
ship dMunari & ZwitteJl 19971) . also noticing that the interstellar 
lines are characterized by multiple components, one of which 
is in emission. In turn, the £V-y=+0.56 value w as calculated 
using the relation given by ICardelli et al.l (119891). hence (V - 
g)o= +0.70 from the 2MASS value g=12.11 (ISkrutskie et al.l 
| 2006l). The latter in dex supplies a T e ff value of 7250 K (Eq. 8 in 
iMasana et alj |2006j). in perfect agreement with the one provided 
by spectroscopy. Table [3] lists the physical parameters obtained 
from space-based photometry and ground-based spectroscopy. 



6. Pulsational model and mode identification 

The parameters obtained from CoRoT photometry and HARPS 
spectroscopy were used to calculate a preliminary pulsational 
model, making certain to match the constraint that f\ must be 
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Table 4. Mode identification from the best pulsational model. 
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Fig. 5. The matching between observed (vertical lines) and calculated 
(points) frequencies for models with degree £ < 3. 



a radial mode. We foll owed the procedur e and the physical as- 
sumptions described in lCasas et al.ld2006l) . In particular, we con- 
sidered a variation in the convective efficiency parameter omlt 
up to 1.5, of the overshoot d m up to 0.3, and an [Fe/H] metallicity 
in the range from -0.52 to 0.08 dex. Other effects such as diffu- 
sion and radiative forces were not included. We assumed a uni- 
form rotation since we could not identify multiplets in advance 
among the 13 frequencies. A non-uniform rotation would be an 
additional source of uncertainty since mixed mod es are highly 
sensitive to th e rotation profile near the stellar core (ISuarez et al.l 
l2006aL[2009h . 

This procedure points out that a further refinement could be 
performed starting from M = 1.71 ±0.02 Mq, L = 18.5+0.2L©, 
R = 2.90 + 0.04%, 7df = 7000 + 50 K, logg = 3.74 + 0.01 dex, 
and [Fe/H]=-0.30+0.20dex. These results indicate a cooler T e ff 
and more pronounced metallic underabundance than those sug- 



gested by the spectroscopic analysis, though the errorbars still 
overlap in both cases. Models with T e g > 7300 K are not able to 
fit f\ as the fundamental radial mode and they moreover iden- 
tify the other frequencies as gravity modes. In a similar way, so- 
lar abundances lead to massive, main-sequence stars. All these 
possibilities are unlikely when considering the whole set of pa- 
rameters and the observed photometric variability. 

We should take into account that the V eq sin i value sug- 
gests a non-negli gible rotational velocity. The adiabatic oscilla- 
tion code FILOU riSuarezfcOOilSuarez & Goupilll2008l and ref- 
erences therein) computes the oscillation frequencies up to the 
second-order effects of rotation. It includes the corrections of 
the frequencies near degeneracy effects occurring when two fre- 
quencies are close to each other. These effects could also occur 
between three or more modes, but they only should be relev ant 
for very high rotational velocities (ISuarez et al.ll2006bll2010h . 

Starting from the reference model pointed out by the spectro- 
scopic and photometric data, we computed sets of models with 
different rotational velocities, ages, and maximum degree t of 
the modes by using the code FILOU. We had to make some as- 
sumptions on the inclination angle i in order to constrain the 
rotational frequency and then infer the rotational splitting. The 
V e g sin i value o f HAPS stars is usually below 25 kins -1 (Fig. 5 
m lBregerll2000l) . suggesting that these variables are intrinsically 
slow rotators. Thus, we considered CoRoT 101155310 as seen 
nearly equator-on, and we allowed variations of V eq in the range 
19-30 km s" 1 . 

The oscillation spectra predicted by the models were com- 
pared with the 1 3 independent terms detected in the frequency 
analysis. The pulsational models limited to I < 3 directly match 
1 1 out of the 13 frequencies. Indeed, modes are predicted around 
all the observed frequencies, except for f$ and fa (Fig.[5j. To fill 
this gap, we computed the models with I < 5 and could obtain 
the mode identifications for these frequencies, too (Table|4]i. The 
average difference between observed and calculated frequencies 
is 0.28 yuHz. The model depends on our assumptions and on the 
uncertainties in the identifications of close frequencies; there- 
fore they cannot be considered unique. We note that identifying 
all frequencies at once with the modes predicted by the models 
with t < 5 improved the fit, because the oscillation spectrum 
would be denser. However, the high prevalence of I > 4 modes 
on £ < 3 ones is contrary to what is expected, since the ampli- 
tudes of the modes should decrease with increasing £ due to the 
cancellation effects. 

The mode identification describes a star very close to the ter- 
minal age on the main sequence, with a central hydrogen content 
X c = 0.15+0.05 and, in turn, an age of 1485 + lOMyr. This evo- 
lutionary stage agrees with the identification of the additional 
terms (i.e., all the terms except f\) with mixed modes. The 
mode is not identified with the first radial overtone, confirming 
that the /1//9 = 0.779 value differs too much from the predicted 
value (0.774). 

We tried a fine tuning of the model after noticing that the 
/3 and /4 terms are identified as two components of a rotation- 
ally split £ = 1 multiplet (Table @J. We considered the existence 
of such multiplet as a further input constraint. To do this, we 
computed models with small variations in the rotational veloc- 
ity around the observed V eq sin i. Such variations had to be large 
enough to modify the m + components sufficiently, and small 
enough to avoid variations in the m = component. The lat- 
ter constraint allows a variation of only about 4 kms~'. Within 
this range, the second-order effects of rotation do not change the 
identification of the other modes. Therefore, the model suggests 
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Fig. 6. Energy balance of the modes of the best fitting model. The green 
box represents where the observed modes are located in frequency. 



that the fz mode definitely does not belong to the same multiplet 
as the f$ and fa modes. 

As a last step, the non-adiabatic analysis of the mode energy 
balance was performed by me ans of the GraCo code dMova et al.l 
l2004HMova & Garridoll2QQ8h . Figure |6] shows the growth range 
of the model. Overstable modes must have a positive growth rate, 
and this happens for a region including the observed frequencies. 

Improvements in the theoretical treatment of the mixed 
modes can lead to a more accurate fitting of the oscillation spec- 
trum of CoRoT 101155310. In any case the mixed modes can 
provide significant constraints on the stellar internal structure 
since now. In particular, our model predicts a convective core 
with size r c ~ 0.08 Rq and a thin outer convective zone lo- 
cated at r - 0.99 jfo, as well as aMLr=0.5, in agreement with 
ICasas et"a fl (l2006h . and d ov = 0.1. These two values affect not 
only the size of the convective zones, but also the possibility of 
estimating the time spent on the main sequence. 



7. Discussion 

The CoRoT satellite offered the possibility of monitoring a vari- 
able star belonging to the class of HADS stars in a continu- 
ous and long way. At the moment, the pulsational behaviour 
disclosed by the timeseries of CoRoT 101155310 is the most 
advanced observational fact we have to understand the physi- 
cal mechanisms at work in the variables of this class. These are 
some questions raised by previous ground-based observations 
on which the CoRoT photometry shed a new light. 

7. 1 . The excitation of independent modes 

The excitation of small-amplitude nonradial modes accompany- 
ing the high-amplitude radial modes in 6 Set variables has been 
debated at length. The pulsational content of V974 Oph provided 
the decisive clear piece of evidence, since this star pulsates in the 
fundamental radial mode (/i=5.23 d _1 , full amplitude 0.32 mag 
in B) and in at least four other independent modes easily de- 
tectable also from ground-based t imeseries (e.g ., /2=5.35 d _1 , 
full amplitude 0. 14 mag in Z?— light: |Porettil2003h . The frequency 
ratios fi/fi of these additional modes (0.734, 0.786, 0.807, and 
0.978) clearly indicate that most of them, if not all, are nonra- 
dial modes. In the case of CoRoT 101 155310 the amplitudes of 
the nonradial modes are much smaller than those of V974 Oph, 
but the superior quality of the CoRoT timeseries does not leave 
any ambiguity in their identification. The ratios f\jf; are listed in 
Table [1] The ratios 0.733, 0.811, and 0.779 are very similar to 
the 0.734, 0.807, and 0.786 ones obtained in the case of V974 



Oph. However, the corresponding modes have a different rank 
in amplitude in the two solutions of the light curves of the two 
stars. 

We can also compare the pulsational content of 
CoRoT 101155310 with that of HD 50844. In this low- 
amplitude 6 Set star, we detected an increase in the signal 
at / < 30 d~' and hundreds of frequencies were neces- 
sary to reduce it to the same level of the noise observed 
at / > 30 d" 1 dPoretfi et al l l2009t) . After removing the 
peaks related to the 61 terms detected in the CoRoT data of 
CoRoT 101155310, the noise level is constant over the full 
5-100 d -1 spectrum (see inserted box in panel e of Fig.0. This 
noise level (29 yU mag) after subtracting 61 frequencies is much 
less than observed in the amplitude spectru m of HD 50844 afte r 
subtracting 500 frequencies (see Fig. 3 in iPoretti et al.l 120091) . 
Moreover, the amplitude of 0. 10 mmag is reached after detecting 
of about 250 frequencies in the light curve of HD 50844, against 
the only 61 ones in the case of CoRoT 101155310. We can 
conclude that the variability of CoRoT 101155310 is much 
simpler than that of HD 50844. New cases will allow us to verify 
if the excitation of only one ten of modes is a characteristic of 
the HADS stars. 



7.2. The modulation of fi 

The physical interpretation of the /„, term is not an easy task. 
We could suppose that it is related with the star's rotation, so 
P=5 . 1 8 1 d is the rotational period. Under this assumption the ra- 
dius value listed in Table[3]supplies an equatorial velocity V eq = 
27 +4 km s" 1 , supporting the hypothesis that CoRoT 101155310 
is seen nearly equator-on. This scenario seems coherent, but we 
have to explain why the rotation affects the amplitude of the main 
pulsation term. A chemical inhomogeneity on the surface could 
do it, but then it should be detected as an independent peak in the 
amplitude spectrum. Moreover, the f m harmonics should also be 
detected, since the equator-on orientation makes the continuous 
visibility unlikely. Actually, frequency f m is not observed in the 
amplitude spectrum of CoRoT 101155310. This lack could be 
due to both its intrinsic, very small amplitude and the instru- 
mental effects that complicate the spectrum at / < 1 d . 

A continuous modulation in amplitude and/or in phase su- 
perposed the main pulsational period is commonly observed in 
RR Lyr stars, and it defines the Blazhko effect. The onset of the 
Blazhko modulation is due to something more specific than ro- 
tation, e.g., the resonance model betw een nonradial, low-degree 
mode s and the main radial mode dDziembowski & Mizerskil 
120041) . or the oblique pulsator model in whi ch the rotational a xis 
does not coincide with the magnetic axis dShiba hashi 2003) or 
the action of a tu rbulent convect i ve dynamo in the lower en- 
velope of the star (IStothersll2006l) . iBregerl (1201 Ol) discusses the 
possibility of explaining the amplitude and phase modulations 
observed in three low-amplitude 5 Set stars (FG Vir, BI CMi, 
and AI CVn) as due to the Blazhko effect. The very rich pulsa- 
tional content of these pulsators favours the interaction between 
modes, without invoking the reson ance mechanism suggested by 
iDziembowski & Mizerskil d2004l) . In the case of HADS stars, 
a much smaller number of modes are excited, and then it ap- 
pears more obvious to link CoRoT 101155310 with RR Lyr stars 
than with low-amplitude 5 Set stars, also having the sharp triplet 
structure in mind. 

It could be useful to investigate the light curves under 
the magnification lens provided by the CoRoT timeseries. The 
long-term stability of the Blazhko modulation in RR Lyr stars 
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is a n ew issue opened u p by the continuous CoRoT time- 
series dPoretti et al.l 120101) . In the case of the HADS variable 
CoRoT 101 155310, there are very low peaks at f\ and harmon- 
ics after subtracting the 61 terms (Fig. |2] panel d), but the in- 
strumental effects revealed by the analysis of the chromatic data 
can be responsible for this residual signal. The amplitudes of f\ 
in white light calculated in each bin are slightly varying during 
the CoRoT monitoring (bottom panel in Fig. 0, but this effect 
is very marginal and could still reflect the complicated changes 
in the trends of the be, gc, and rc data. The light curves of 
Blazhko RR Lyr stars show b umps and larger sc atter at well- 
defined phases (e.g., Fig. 3 in iPoretti et al.ll20loh . probably in 
relation with shock waves. In the case of CoRoT 101155310, 
the smoothed light curve over /j (Fig. [3] upper curve in the top 
panel) and the residuals folded with f\ (Fig. [3] lower curve in the 
top panel) do not show any particularity. Therefore, at the mo- 
ment we have no further evidence supporting the straight iden- 
tification of the triplets around f\ as due to the Blazhko effect 
observed in RR Lyr stars. 

In RR Lyr stars, the amplitude of the Blazhko modulation is 
stronger than that of additional modes. Indeed, CoRoT photom- 
etry allowed us to study the Blazhko effect in great detail and the 
discovery of addit ional modes was a bonus dChadid et al.ll2010t 
IPoretti et al.ll2010l) . In the case of HADS stars, the opposite case 
occurred; i.e., searching for additional modes we discovered the 
fx modulation. This demonstrates once more how rewarding the 
investment in the observation of known phenomena with more 
advanced tools could be. 

8. Conclusions 

The solution for the CoRoT light curve of CoRoT 101155310 
allowed us to detect only 12 additional modes; i.e., we did not 
detect the huge number of fr equencies found in low-amplitude 
pulsators dPoretti et al.ll2009b . These modes all have mixed na- 
ture and, together with the fundamental radial mode, originate 
sixty-one frequencies able to solve the CoRoT light curve. We 
can speculate that the large amount of energy put in the radial 
modes limits the number of excited nonradial modes, but the 
observations of new HADS stars could shed new light on this 
particular aspect. 

The newest result is the small but very clear modulation of 
the fundamental radial mode, observed for the first time in an 
HADS star. The observation of other HADS stars by space pho- 
tometry could tell us whether this is a rotational effect or an 
extension of the Blazhko effect from RR Lyr to HADS stars. 
The similarity between these two classes of variables would be a 
new way of looking at the Hertzsprung-Russell diagram and, if 
confirmed, would be able to put tighter constraints on the phys- 
ical model of the Blazhko effect, also taking the thin convec- 
tive layer predicted by the model of CoRoT 101 155310 into ac- 
count. Indeed, it should be able to explain the modulations ob- 
served both in low-mass stars in a core He-burning phase (RR 
Lyr stars) and in more massive, post-main sequence stars moving 
to the giant branch (HADS stars). In this context, nonradial res- 
onant modes close to the radial fundamental mode seem to have 
the requisites needed to act in such different evolutionary stages, 
but the other models are still plausible. We also note that the 
modulation observed in CoRoT 101155310 (5.181 d) is shorter 
in time than the sho rtest ever observed for Blazhko RR Lyr stars 
(5.309 d in SS Cnc: lJurcsik e"t~aT1l2006l) . 

The excitation of mixed modes in stars approaching the ter- 
minal age on the main sequence is confirmed and they are found 
in the predicted overstable region. Further improvements in the 



theoretical modelling are necessary, especially when dealing 
with close frequencies in the presence of a non-negligible ro- 
tation. 
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